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Abstract —An analytic model is proposed of a solar power plant (SPP) with a Stirling engine that is based on 
the isothermal model of the Stirling engine (SE) working process and is improved by account for the actual heat 
exchange, the hydraulic and mechanical losses in the SE, the losses in the electric generator, and also the basic 
parameters of the solar radiation concentrator. 
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Theoretically, the ideal Stirling cycle has the maxi¬ 
mal possible efficiency for a given temperature differ¬ 
ential, equal to the efficiency of the Carnot cycle. How¬ 
ever the real working process in the Stirling engine dif¬ 
fers from the ideal process, and its effectiveness is 
considerably lower. This is due to the irreversibility of 
the actual thermodynamic processes, the deviation of 
the working processes in the working chambers of the 
Stirling engine from isothermality, the nonideality of 
the heat exchange processes in the heater, the regenera¬ 
tor, and the cooler, and also the heat losses to the ambi¬ 
ent medium. 

The program for calculating the characteristics of 
the Stirling engine is based on the Schmidt isothermal 
model that is improved by accounting for the heat trans¬ 
fer, the hydraulic resistances, the nonideality of the 
regenerator, the mechanical losses, and the losses in the 
electric generator [1, 2], Moreover, the sinusoidal 
motion of the pistons is replaced by the equations of the 
actual motion of the pistons with account for the spe¬ 
cific drive mechanism. In the calculation of the solar 
power plant as a whole, the losses in the concentrator- 
receiver system are taken into account. 

Figure 1 presents the scheme of the 5-kW Stirling 
engine with a swashplate drive that was developed at 
FTI NPO “Fizika-Solntse” of the Academy of Sciences, 
Republic of Uzbekistan. The algorithm for the calcula¬ 
tion of the engines of this type is presented later. 

The piston travel is defined by the expression 

h = 2Rssin(a), (1) 

where a is the swashplate inclination angle; and Rs is 
the radius of the swashplate. 

The instantaneous piston position 1 is defined by the 
expression 


The reference point is to t = 0, the piston is at the 
BDC point, and VG = max. 

The volume of the hot chamber of the first cylinder is 

Vg = 0.785 Dcl 2 (Z i + Hgm)\ (3) 

the volume of the cold chamber of the adjacent (sec¬ 
ond) cylinder is 

V, = 0.785 (Del 2 - Dst 2 )(h + Hxm - Z 2 ), (4) 

where Hgm and Hxm are the piston clearances at the 
TDC and BDC points; Del is the diameter of the cylin¬ 
der; Dst is the diameter of the rod; and Z 2 = 0.5/i(l + 
cos(cof + 7t/2)). 

The calculation of the working fluid pressure func¬ 
tion for the given initial temperatures of the working 



Z x - 0.5/i(l + cos(cof)). 


( 2 ) 


Fig. 1 . Solar Stirling engine with swashplate drive. 
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chambers with the absence of hydraulic resistances in tudes of the instantaneous and cycle-average heat trans- 
the heat exchangers is performed as follows: fer coefficients in the heater and the cooler: 


P((Ot) = Psr(l -A')/(l + Acos(cot-fi)), (5) 

where A x - (x 2 + K 2 + 2x7fcos(7t/2)) 0 - 5 ; A 2 = X + K + 
2 (Xln + Xc + 2(Xn)/(\ + x)); A = A,/A 2 ; B { = 
Ksir\(n/2)/AB = Asin(P,); x = Txs/Tgs; Tgs, Txs are 
the temperatures of the walls of the heater and the 
cooler. 

The calculation of the working fluid masses in the 
hot and cold chambers and in the heat exchangers of the 
engine is performed on the basis of the condition that 
the working fluid is an ideal gas: 

Mg(cot) = P( cot) Vg(cot)/( RTg) is the working fluid 
mass in the hot chamber; 

Mli(o'>t)P( (m) V| Xh/(RTg) is the working fluid mass 
in the heater; 

Mr((Ot) = 2P((Ot)V 1 Xr/(R(Tg + Tx)) is the working 
fluid mass in the regenerator; 

Mc(cot)P(cot) V t Xc/(RTx) is the working fluid mass 
in the cooler; 

Mx(cot) = P( OK) Vx( I)/(RTx) is the working fluid 
mass in the cold chamber. 

The calculation of the heat transfer coefficients in 
the heat exchangers (heater and cooler) is based on the 
assumption of quasistationarity of the heat exchange 
processes, for the instantaneous values of the gas veloc¬ 
ity in the heat exchanger channels, 


a,,(cot) = Nuh{(f)t)XJDIr, 

( 12 ) 

oq(CGt) = Nuc(tot)X ( ./Z)c; 
a hsr = Y a/,(cot)/(27t/Acot); 

(13) 

a csr = ^a £ .(cot)/(27t/Acot). 

For the calculation of the average temperatures in 
the working chambers of the engine, we arbitrarily 
specify the initial magnitude of the heat that is supplied 
to the cycle: 

2, = 0A5PsrNV x . (14) 

Then the magnitudes of the temperature differen¬ 
tials on the heat exchange surfaces arc defined by the 
expressions 

ATh = Q x Dhl{a h AXhV x )\ (15) 


A 7c = Q x xDd(a c 4XcV x ). (16) 

The gas temperatures in the hot and cold chambers 
of the engine are determined by the relations Tg = Tg.v - 
AT/?, Tx = Tx.v - AT/?, and their ratio = Tx/Tg. 

The working fluid pressure functions in the cold and 
hot chambers of the engine with account for the real 
heat exchange and the actual hydraulic resistances are 
determined by the expressions 


Reh(cot) = AcotNLhDhABS(dMg(cot) 

+ 0.5dMh(cot))/(V l [ii,)-, 

(6) 

Px(oot) = 7 > (cot)(l - AP(cot) 

X ((1 + cos(C 0 t))xr + 2Xhlr 

(17) 

Rec(cot) = AcotNLcDcABS(dMx(cot ) 

+ 0.5dMc(I))/(V ] Hc)- 

(7) 

+ 2XrXr/( 1 + xr))/(A 2 P.vr)); 

Pg(cot) = Px(cot) + AP(cot), 

(18) 


For the laminar flow regime with Re < 2200, the 
Nusselt number for each heat exchanger is calculated 
using the expression [3] 

Nu(cot) = 0.142Re(cot) 033 . (8) 

For the transition interval, the formulas for calculat¬ 
ing the Nusselt number are derived from the tabulated 
values: 

for Re(cot) > 2200 and < 5600 

Nu(cot) = 32.8 x 0.43429ALOG 

(9) 

x (Re(cof)) - 107.57; 

for Re(cot) > 5600 and < 10000 

Nu(cot) = 48.07 x 0.43429ALOG 

( 10 ) 

x(Re(cot))- 164.88. 

If Re(cot) > 10000, i.e., for the turbulent regime 

Nu(cot) = 0.0175(Re(cot))° 8 . (11) 

The calculation of the heat transfer in the cooler is 
performed similarly. As a result, we obtain the magni- 


where AP(cot) is the overall hydraulic resistance in the 
channels of the heat exchangers and the connecting 
lines, which can be determined from the known rela¬ 
tions for the instantaneous working fluid velocities. 

The calculation of the indicated work in the hot and 
cold chambers of the engine: 

Lig = J^(®t)dV'g(cot); (19) 

Lix = jPx((Ot)dVx((Ot). (20) 

The regenerator effectiveness T|r is determined 
using the technique described in [4]. 

The heat that is supplied to the hot chamber of the 
engine is determined by the expression 

Q = Lig+ Qr{ \ -T|r), (21) 

where Qr = CpGr(Tg - Tx) is the heat that is supplied 
to and removed from the regenerator per cycle, and 
Gr is the gas flowrate in the regenerator. 

Then the calculation is performed by the successive 
approximation method. For example, if ABS((9| - Q)Q\ > 
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> 0.001, i.e., the initial and obtained values differ by a 
magnitude that is larger than 0.001, the calculation is 
repeated until the accuracy of the calculation satisfies 
the specified condition. In this case the value (Q x + Q)/2 
is specified as the initial value Q x . As a result, we 
obtain: 

the indicated work per cycle Li = Lig + Lix; (22) 
the indicated power Ni = Lin ; (23) 

the indicated engine efficiency r\i - LHQ. (24) 

For the presentation of the results of the calcula¬ 
tions, it is helpful to use the specific power, which is 
defined by the expression 

PM = M/(PSR(1 + A)V». (25) 

We shall use the dimensionless power for the calcu¬ 
lation of the indicated power and the determination of 
the parameters of the heat exchangers for the engines of 
analogous design, in the present case of the alpha type, 
but with differing working volume. 

The friction process in the friction assemblies of the 
cylinder-piston group of the Stirling engine takes place 
with reciprocating motion of the piston rings, which are 
made of antifriction materials, relative to the steel sur¬ 
face of the cylinder in the conditions of the absence of 
a lubricant. A combination of guide and sealing rings is 
usually installed on the pistons. The friction forces in 
the sealing rings of the SE arise because they are 
pushed against the cylinder wall by the forces of the gas 
pressure and by the elastic force of the deformation of 
the expander in the rectilinear reciprocating motion of 
the piston. The friction forces in the guide rings arise 
because they are pushed against the cylinder wall by the 
lateral components of the forces from the drive mecha¬ 
nism. 

The mechanical losses on friction in the piston rings 
are calculated using the technique described in [5], 

The mechanical power on the engine shaft with 
account for the losses in the piston-cylinder group and 
the drive is defined by the expression 

Ne = 4(M-ATr)r| dr . (26) 

The mechanical (effective) efficiency of the engine 
is determined as 

r\e = Ne/(4 Q). (27) 

The electric power and efficiency of the engine are: 

Nel = Ne/r| gen ; (28) 

pel = Nel/(4<2), (29) 

where p gen =/(Nel). The study of the characteristics of 
the electric generator was performed earlier on the one- 
kW and five-kW power levels [6, 7], 

The efficiency of the solar power plant as a whole is 
defined by the ratio of the generated electric power to 
the solar flux incident on the concentrator: 


0 pad , Q pr , Nel, kW 



Fig. 2. Dependences of the power of the solar radiation flux 
that is incident on the concentrator, the thermal power that 
is transferred to the engine working fluid, and the generated 
electric power on the density of the direct solar radiation. 


T| spp = Nel/(EsFk). (30) 

For the calculation of the required concentrator area, 
it is necessary to consider the losses due to the accuracy 
of the concentrator tracking system, the alignment 
losses, the heat receiver losses, and also the density of 
the direct solar radiation. 

The calculation of the torque is performed for each 
cylinder, and then, summing these values for each 
moment of time, we obtain the torque on the engine 
shaft. The force that acts on the shaft from a single pis¬ 
ton is 

Fkri(cot) = AP/(cot)tn(a)(-sin(cot)) 5 (31) 

and the torque from a single piston is 
Mkri(cot) 

Fkri(cot)Rscos(a) (32) 

Jsin(a)t + 0.57t) 2 + cos(a)“cos(cot + 0.57t) 2 

The torque on the engine shaft is determined by the 
sum of the instantaneous values from all the pistons: 

Mkr(cot) = ^ Mkri(cot ). (33) 

The calculations were made for a solar power plant 
with a 5-kW Stirling engine. The following characteris¬ 
tics of the concentrator were adopted for the calcula¬ 
tion: the concentrator area is 40 m 2 ; the reflection coef¬ 
ficient is 0.85; the alignment losses and the tracking 
inaccuracy losses are 5%. The results of the calculation 
are presented in Figs. 2 and 3. 

The results of the calculation agree with the experi¬ 
mental data that were obtained in the tests of a solar 
power plant with a V-160 engine rated at 7.5 kW. The 
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SPP efficiency 



Fig. 3. Dependence of efficiency of SPP with 5-kW Stirling 
engine on the density of the direct solar radiation. 


tests were performed at the solar test station at Almeria 
(Spain) [8]. The developed analytic technique can be used 
for the analysis and preparation of the basic data in the 
design of the solar power plant with a Stirling engine. 
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